A pan-cancer clinical platform to predict immunotherapy outcomes and prioritize immuno-oncology combinations in early-phase trials.

BACKGROUND: Immunotherapy is effective, but current biomarkers for patient selection have proven modest sensitivity. Here, we developed VIGex, an optimized gene signature based on the expression level of 12 genes involved in immune response with RNA sequencing. METHODS: We implemented VIGex using the nCounter platform (Nanostring) on a large clinical cohort encompassing 909 tumor samples across 45 tumor types. VIGex was developed as a continuous variable, with cutoffs selected to detect three main categories (hot, intermediate-cold and cold) based on the different inflammatory status of the tumor microenvironment. FINDINGS: Hot tumors had the highest VIGex scores and exhibited an increased abundance of tumor-infiltrating lymphocytes as compared with the intermediate-cold and cold. VIGex scores varied depending on tumor origin and anatomic site of metastases, with liver metastases showing an immunosuppressive tumor microenvironment. The predictive power of VIGex-Hot was observed in a cohort of 98 refractory solid tumor from patients treated in early-phase immunotherapy trials and its clinical performance was confirmed through an extensive metanalysis across 13 clinically annotated gene expression datasets from 877 patients treated with immunotherapy agents. Last, we generated a pan-cancer biomarker platform that integrates VIGex categories with the expression levels of immunotherapy targets under development in early-phase clinical trials. CONCLUSIONS: Our results support the clinical utility of VIGex as a tool to aid clinicians for patient selection and personalized immunotherapy interventions. FUNDING: BBVA Foundation; 202-2021 Division of Medical Oncology and Hematology Fellowship award; Princess Margaret Cancer Center.

Human Metastatic Cholangiocarcinoma Patient-Derived Xenografts and Tumoroids for Preclinical Drug Evaluation.

PURPOSE: Cholangiocarcinoma (CCA) is usually diagnosed at advanced stages, with limited therapeutic options. Preclinical models focused on unresectable metastatic CCA are necessary to develop rational treatments. Pathogenic mutations in IDH1/2, ARID1A/B, BAP1, and BRCA1/2 have been identified in 30%-50% of patients with CCA. Several types of tumor cells harboring these mutations exhibit homologous recombination deficiency (HRD) phenotype with enhanced sensitivity to PARP inhibitors (PARPi). However, PARPi treatment has not yet been tested for effectiveness in patient-derived models of advanced CCA. EXPERIMENTAL DESIGN: We have established a collection of patient-derived xenografts from patients with unresectable metastatic CCA (CCA_PDX). The CCA_PDXs were characterized at both histopathologic and genomic levels. We optimized a protocol to generate CCA tumoroids from CCA_PDXs. We tested the effects of PARPis in both CCA tumoroids and CCA_PDXs. Finally, we used the RAD51 assay to evaluate the HRD status of CCA tissues. RESULTS: This collection of CCA_PDXs recapitulates the histopathologic and molecular features of their original tumors. PARPi treatments inhibited the growth of CCA tumoroids and CCA_PDXs with pathogenic mutations of BRCA2, but not those with mutations of IDH1, ARID1A, or BAP1. In line with these findings, only CCA_PDX and CCA patient biopsy samples with mutations of BRCA2 showed RAD51 scores compatible with HRD. CONCLUSIONS: Our results suggest that patients with advanced CCA with pathogenic mutations of BRCA2, but not those with mutations of IDH1, ARID1A, or BAP1, are likely to benefit from PARPi therapy. This collection of CCA_PDXs provides new opportunities for evaluating drug response and prioritizing clinical trials.

The trophectoderm acts as a niche for the inner cell mass through C/EBPα-regulated IL-6 signaling.

IL-6 has been shown to be required for somatic cell reprogramming into induced pluripotent stem cells (iPSCs). However, how Il6 expression is regulated and whether it plays a role during embryo development remains unknown. Here, we describe that IL-6 is necessary for C/EBPα-enhanced reprogramming of B cells into iPSCs but not for B cell to macrophage transdifferentiation. C/EBPα overexpression activates both Il6 and Il6ra genes in B cells and in PSCs. In embryo development, Cebpa is enriched in the trophectoderm of blastocysts together with Il6, while Il6ra is mostly expressed in the inner cell mass (ICM). In addition, Il6 expression in blastocysts requires Cebpa. Blastocysts secrete IL-6 and neutralization of the cytokine delays the morula to blastocyst transition. The observed requirement of C/EBPα-regulated IL-6 signaling for pluripotency during somatic cell reprogramming thus recapitulates a physiologic mechanism in which the trophectoderm acts as niche for the ICM through the secretion of IL-6.

The HDAC7-TET2 epigenetic axis is essential during early B lymphocyte development.

Correct B cell identity at each stage of cellular differentiation during B lymphocyte development is critically dependent on a tightly controlled epigenomic landscape. We previously identified HDAC7 as an essential regulator of early B cell development and its absence leads to a drastic block at the pro-B to pre-B cell transition. More recently, we demonstrated that HDAC7 loss in pro-B-ALL in infants associates with a worse prognosis. Here we delineate the molecular mechanisms by which HDAC7 modulates early B cell development. We find that HDAC7 deficiency drives global chromatin de-condensation, histone marks deposition and deregulates other epigenetic regulators and mobile elements. Specifically, the absence of HDAC7 induces TET2 expression, which promotes DNA 5-hydroxymethylation and chromatin de-condensation. HDAC7 deficiency also results in the aberrant expression of microRNAs and LINE-1 transposable elements. These findings shed light on the mechanisms by which HDAC7 loss or misregulation may lead to B cell-based hematological malignancies.

Coding and non-coding co-expression network analysis identifies key modules and driver genes associated with precursor lesions of gastric cancer.

BACKGROUND: Helicobacter pylori infection is the most important risk factor for gastric cancer (GC). Human gastric adenocarcinoma develops after long-term H. pylori infection via the Correa cascade. This carcinogenic pathway describes the progression from gastritis to atrophy, intestinal metaplasia (IM), dysplasia and GC. Patients with atrophy and intestinal metaplasia are considered to have precancerous lesions of GC (PLGC). H. pylori eradication and endoscopy surveillance are currently the only interventions for preventing GC. Better knowledge of the biology of human PLGC may help find stratification markers and contribute to better understanding of biological mechanisms. One way to achieve this is by using co-expression network analysis. Weighted gene co-expression network analysis (WGCNA) is often used to identify modules from co-expression networks and relate them to clinical traits. It also allows identification of driver genes that may be critical for PLGC. AIM: The purpose of this study was to identify co-expression modules and differential gene expression in dyspeptic patients at different stages of the Correa pathway. METHODS: We studied 96 gastric biopsies from 78 patients that were clinically classified as: non-active (n = 10) and chronic-active gastritis (n = 20), atrophy (n = 12), and IM (n = 36). Gene expression of coding RNAs was determined by microarrays and non-coding RNAs by RNA-seq. The WGCNA package was used for network construction, module detection, module preservation and hub and driver gene selection. RESULTS: WGCNA identified 20 modules for coding RNAs and 4 for each miRNA and small RNA class. Modules were associated with antrum and corpus gastric locations, chronic gastritis and IM. Notably, coding RNA modules correlated with the Correa cascade. One was associated with the presence of H. pylori. In three modules, the module eigengene (ME) gradually increased in the stages toward IM, while in three others the inverse relationship was found. One miRNA module was negatively correlated to IM and was used for a mRNA-miRNA integration analysis. WGCNA also uncovered driver genes. Driver genes show both high connectivity within a module and are significantly associated with clinical traits. Some of those genes have been previously involved in H. pylori carcinogenesis, but others are new. Lastly, using similar external transcriptomic data, we confirmed that the discovered mRNA modules were highly preserved. CONCLUSION: Our analysis captured co-expression modules that provide valuable information to understand the pathogenesis of the progression of PLGC.

MicroRNAs Deregulated in Intraductal Papillary Mucinous Neoplasm Converge on Actin Cytoskeleton-Related Pathways That Are Maintained in Pancreatic Ductal Adenocarcinoma.

Intraductal papillary mucinous neoplasms (IPMN) are pancreatic cystic lesions that can develop into pancreatic ductal adenocarcinoma (PDAC). Although there is an increasing incidence of IPMN diagnosis, the mechanisms of formation and progression into invasive cancer remain unclear. MicroRNAs (miRNAs) are small non-coding RNAs, repressors of mRNA translation, and promising diagnostic biomarkers for IPMN and PDAC. Functional information on the role of early-altered miRNAs in this setting would offer novel strategies for tracking the IPMN-to-PDAC progression. In order to detect mRNAs that are likely to be under miRNA regulation in IPMNs, whole transcriptome and miRNome data from normal pancreatic tissue (n = 3) and IPMN lesions (n = 4) were combined and filtered according to negative correlation and miRNA-target prediction databases by using miRComb R package. Further comparison analysis with PDAC data allowed us to obtain a subset of miRNA-mRNA pairs shared in IPMN and PDAC. Functional enrichment analysis unravelled processes that are mainly related with cell structure, actin cytoskeleton, and metabolism. MiR-181a appeared as a master regulator of these processes. The expression of selected miRNA-mRNA pairs was validated by qRT-PCR in an independent cohort of patients (n = 40), and then analysed in different pancreatic cell lines. Finally, we generated a cellular model of HPDE cells stably overexpressing miR-181a, which showed a significant alteration of actin cytoskeleton structures accompanied by a significant downregulation of EPB41L4B and SEL1L expression. In situ hybridization of miR-181a and immunohistochemistry of EPB41L4B and SEL1L in pancreatic tissues (n = 4 Healthy; n = 3 IPMN; n = 4 PDAC) were also carried out. In this study, we offer insights on the potential implication of miRNA alteration in the regulation of structural and metabolic changes that pancreatic cells experience during IPMN establishment and that are maintained in PDAC.

Colon capsule endoscopy versus CT colonography in FIT-positive colorectal cancer screening subjects: a prospective randomised trial-the VICOCA study.

BACKGROUND: Colon capsule endoscopy (CCE) and CT colonography (CTC) are minimally invasive techniques for colorectal cancer (CRC) screening. Our objective is to compare CCE and CTC for the identification of patients with colorectal neoplasia among participants in a CRC screening programme with positive faecal immunochemical test (FIT). Primary outcome was to compare the performance of CCE and CTC in detecting patients with neoplastic lesions. METHODS: The VICOCA study is a prospective, single-centre, randomised trial conducted from March 2014 to May 2016; 662 individuals were invited and 349 were randomised to CCE or CTC before colonoscopy. Endoscopists were blinded to the results of CCE and CTC. RESULTS: Three hundred forty-nine individuals were included: 173 in the CCE group and 176 in the CTC group. Two hundred ninety individuals agreed to participate: 147 in the CCE group and 143 in the CTC group. In the intention-to-screen analysis, sensitivity, specificity and positive and negative predictive values for the identification of individuals with colorectal neoplasia were 98.1%, 76.6%, 93.7% and 92.0% in the CCE group and 64.9%, 95.7%, 96.8% and 57.7% in the CTC group. In terms of detecting significant neoplastic lesions, the sensitivity of CCE and CTC was 96.1% and 79.3%, respectively. Detection rate for advanced colorectal neoplasm was higher in the CCE group than in the CTC group (100% and 93.1%, respectively; RR = 1.07; p = 0.08). Both CCE and CTC identified all patients with cancer. CCE detected more patients with any lesion than CTC (98.6% and 81.0%, respectively; RR = 1.22; p = 0.002). CONCLUSION: Although both techniques seem to be similar in detecting patients with advanced colorectal neoplasms, CCE is more sensitive for the detection of any neoplastic lesion. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT02081742 . Registered: September 16, 2013.

Cross-sectional study of human coding- and non-coding RNAs in progressive stages of Helicobacter pylori infection.

Helicobacter pylori infects 4.4 billion individuals worldwide and is considered the most important etiologic agent for peptic ulcers and gastric cancer. Individual response to H. pylori infection is complex and depends on complex interactions between host and environmental factors. The pathway towards gastric cancer is a sequence of events known as Correa's model of gastric carcinogenesis, a stepwise inflammatory process from normal mucosa to chronic-active gastritis, atrophy, metaplasia and gastric adenocarcinoma. This study examines gastric clinical specimens representing different steps of the Correa pathway with the aim of identifying the expression profiles of coding- and non-coding RNAs that may have a role in Correa's model of gastric carcinogenesis. We screened for differentially expressed genes in gastric biopsies by employing RNAseq, microarrays and qRT-PCR. Here we provide a detailed description of the experiments, methods and results generated. The datasets may help other scientists and clinicians to find new clues to the pathogenesis of H. pylori and the mechanisms of progression of the infection to more severe gastric diseases. Data is available via ArrayExpress.

Copy-Number Alteration Burden Differentially Impacts Immune Profiles and Molecular Features of Hepatocellular Carcinoma.

PURPOSE: Chromosomal instability is a hallmark of cancer that results in broad and focal copy-number alterations (CNAs), two events associated with distinct molecular, immunologic, and clinical features. In hepatocellular carcinoma (HCC), the role of CNAs has not been thoroughly assessed. Thus, we dissected the impact of CNA burdens on HCC molecular and immune features. EXPERIMENTAL DESIGN: We analyzed SNP array data from 452 paired tumor/adjacent resected HCCs and 25 dysplastic nodules. For each sample, broad and focal CNA burdens were quantified using CNApp, and the resulting broad scores (BS) and focal scores (FS) were correlated with transcriptomic, mutational, and methylation profiles, tumor immune composition, and clinicopathologic data. RESULTS: HCCs with low broad CNA burdens (defined as BS ≤ 4; 17%) presented high inflammation, active infiltrate signaling, high cytolytic activity, and enrichment of the "HCC immune class" and gene signatures related to antigen presentation. Conversely, tumors with chromosomal instability (high broad CNA loads, BS ≥ 11; 40%), displayed immune-excluded traits and were linked to proliferation, TP53 dysfunction, and DNA repair. Candidate determinants of the low cytotoxicity and immune exclusion features of high-BS tumors included alterations in antigen-presenting machinery (i.e., HLA), widespread hypomethylation, and decreased rates of observed/expected neoantigenic mutations. High FSs were independent of tumor immune features, but were related to proliferation, TP53 dysfunction, and progenitor cell traits. CONCLUSIONS: HCCs with high chromosomal instability exhibit features of immune exclusion, whereas tumors displaying low burdens of broad CNAs present an immune active profile. These CNA scores can be tested to predict response to immunotherapies.

CTCF is dispensable for immune cell transdifferentiation but facilitates an acute inflammatory response.

Three-dimensional organization of the genome is important for transcriptional regulation1-7. In mammals, CTCF and the cohesin complex create submegabase structures with elevated internal chromatin contact frequencies, called topologically associating domains (TADs)8-12. Although TADs can contribute to transcriptional regulation, ablation of TAD organization by disrupting CTCF or the cohesin complex causes modest gene expression changes13-16. In contrast, CTCF is required for cell cycle regulation17, embryonic development and formation of various adult cell types18. To uncouple the role of CTCF in cell-state transitions and cell proliferation, we studied the effect of CTCF depletion during the conversion of human leukemic B cells into macrophages with minimal cell division. CTCF depletion disrupts TAD organization but not cell transdifferentiation. In contrast, CTCF depletion in induced macrophages impairs the full-blown upregulation of inflammatory genes after exposure to endotoxin. Our results demonstrate that CTCF-dependent genome topology is not strictly required for a functional cell-fate conversion but facilitates a rapid and efficient response to an external stimulus.

MiR-93 is related to poor prognosis in pancreatic cancer and promotes tumor progression by targeting microtubule dynamics.

Biomarkers and effective therapeutic agents to improve the dismal prognosis of pancreatic ductal adenocarcinoma (PDAC) are urgently required. We aimed to analyze the prognostic value and mechanistic action of miR-93 in PDAC. Correlation of miR-93 tumor levels from 83 PDAC patients and overall survival (OS) was analyzed by Kaplan-Meier. MiR-93 depletion in PANC-1 and MIA PaCa-2 cells was achieved by CRISPR/Cas9 and miR-93 overexpression in HPDE cells by retroviral transduction. Cell proliferation, migration and invasion, cell cycle analysis, and in vivo tumor xenografts in nude mice were assessed. Proteomic analysis by mass spectrometry and western-blot was also performed. Finally, miR-93 direct binding to candidate mRNA targets was evaluated by luciferase reporter assays. High miR-93 tumor levels are significantly correlated with a worst prognosis in PDAC patients. MiR-93 abolition altered pancreatic cancer cells phenotype inducing a significant increase in cell size and a significant decrease in cell invasion and proliferation accompanied by a G2/M arrest. In vivo, lack of miR-93 significantly impaired xenograft tumor growth. Conversely, miR-93 overexpression induced a pro-tumorigenic behavior by significantly increasing cell proliferation, migration, and invasion. Proteomic analysis unveiled a large group of deregulated proteins, mainly related to G2/M phase, microtubule dynamics, and cytoskeletal remodeling. CRMP2, MAPRE1, and YES1 were confirmed as direct targets of miR-93. MiR-93 exerts oncogenic functions by targeting multiple genes involved in microtubule dynamics at different levels, thus affecting the normal cell division rate. MiR-93 or its direct targets (CRMP2, MAPRE1, or YES1) are new potential therapeutic targets for PDAC.

CNApp, a tool for the quantification of copy number alterations and integrative analysis revealing clinical implications.

Somatic copy number alterations (CNAs) are a hallmark of cancer, but their role in tumorigenesis and clinical relevance remain largely unclear. Here, we developed CNApp, a web-based tool that allows a comprehensive exploration of CNAs by using purity-corrected segmented data from multiple genomic platforms. CNApp generates genome-wide profiles, computes CNA scores for broad, focal and global CNA burdens, and uses machine learning-based predictions to classify samples. We applied CNApp to the TCGA pan-cancer dataset of 10,635 genomes showing that CNAs classify cancer types according to their tissue-of-origin, and that each cancer type shows specific ranges of broad and focal CNA scores. Moreover, CNApp reproduces recurrent CNAs in hepatocellular carcinoma and predicts colon cancer molecular subtypes and microsatellite instability based on broad CNA scores and discrete genomic imbalances. In summary, CNApp facilitates CNA-driven research by providing a unique framework to identify relevant clinical implications. CNApp is hosted at https://tools.idibaps.org/CNApp/.

In most cases, human cells contain two copies of each of their genes, yet sometimes this can change, an effect called copy number alteration (CNA). Cancer is a genetic disease and thus, studying the DNA from tumor samples is crucial to improving diagnosis and choosing the right treatment. Most tumors contain cells with CNAs; however, the impact of CNAs in cancer progression is poorly understood. CNAs can be studied by examining the genome of tumor cells and finding which regions display an unusual number of copies. It may also be possible to gather information about different cancer types by analyzing the CNAs in a tumor, but this approach requires the analysis of large amounts of data. To aid the analysis of CNAs in cancer cells, Franch-Expósito, Bassaganyas et al. have created an online tool called CNApp, which is able to identify and count CNAs in genomic data and link them to features associated with different cancers. The hope is that a better understanding of the effect of CNAs in cancer could help better diagnose cancers, and improve outcomes for patients. Potentially, this could also predict what type of treatment would work better for a specific tumor. Besides, by using a machine-learning approach, the tool can also make predictions about specific cancer subtypes in order to facilitate clinical decisions. Franch-Expósito, Bassaganyas et al. tested CNApp using previously existing cancer data from 33 different cancer types to show how CNApp can help the interpretation of CNAs in cancer. Moreover, CNApp can also use CNAs to identify different types of bowel (colorectal) cancer in a way that could help doctors to make decisions about treatment. Together these findings show that CNApp provides an adaptable and accessible research tool for the study of cancer genomics, which could provide opportunities to inform medical procedures.

Whsc1 links pluripotency exit with mesendoderm specification.

How pluripotent stem cells differentiate into the main germ layers is a key question of developmental biology. Here, we show that the chromatin-related factor Whsc1 (also known as Nsd2 and MMSET) has a dual role in pluripotency exit and germ layer specification of embryonic stem cells. On induction of differentiation, a proportion of Whsc1-depleted embryonic stem cells remain entrapped in a pluripotent state and fail to form mesendoderm, although they are still capable of generating neuroectoderm. These functions of Whsc1 are independent of its methyltransferase activity. Whsc1 binds to enhancers of the mesendodermal regulators Gata4, T (Brachyury), Gata6 and Foxa2, together with Brd4, and activates the expression of these genes. Depleting each of these regulators also delays pluripotency exit, suggesting that they mediate the effects observed with Whsc1. Our data indicate that Whsc1 links silencing of the pluripotency regulatory network with activation of mesendoderm lineages.

Induced Chromosomal Aneuploidy Results in Global and Consistent Deregulation of the Transcriptome of Cancer Cells.

Chromosomal aneuploidy is a defining feature of epithelial cancers. The pattern of aneuploidies is cancer-type specific. For instance, the gain of chromosome 13 occurs almost exclusively in colorectal cancer. We used microcell-mediated chromosome transfer to generate gains of chromosome 13 in the diploid human colorectal cancer cell line DLD-1. Extra copies of chromosome 13 resulted in a significant and reproducible up-regulation of transcript levels of genes on chromosome 13 (P = .0004, FDR = 0.01) and a genome-wide transcriptional deregulation in all 8 independent clones generated. Genes contained in two clusters were particularly affected: the first cluster on cytoband 13q13 contained 7 highly up-regulated genes (NBEA, MAB21L1, DCLK1, SOHLH2, CCDC169, SPG20 and CCNA1, P = .0003) in all clones. A second cluster was located on 13q32.1 and contained five upregulated genes (ABCC4, CLDN10, DZIP1, DNAJC3 and UGGT2, P = .003). One gene, RASL11A, localized on chromosome band 13q12.2, escaped the copy number-induced overexpression and was reproducibly and significantly down-regulated on the mRNA and protein level (P = .0001, FDR = 0.002). RASL11A expression levels were also lower in primary colorectal tumors as compared to matched normal mucosa (P = .0001, FDR = 0.0001. Overexpression of RASL11A increases cell proliferation and anchorage independent growth while decreasing cell migration in +13 clones. In summary, we observed a strict correlation of genomic copy number and resident gene expression levels, and aneuploidy dependent consistent genome-wide transcriptional deregulation.

Novel Circulating miRNA Signatures for Early Detection of Pancreatic Neoplasia.

OBJECTIVES: Pancreatic ductal adenocarcinoma (PDAC) presents the lowest survival rate of all cancers because only 6% of patients reach five-year survival. Alterations in the expression of several microRNAs (miRNAs) occur in the tumor of PDAC and in preneoplastic lesions as the called intraductal papillary mucinous neoplasm (IPMN). Here, we aimed at identifying which miRNAs are significantly altered in liquid biopsies from patients with PDAC and IPMN to find new noninvasive biomarkers for early detection of PDAC. METHODS: We analyzed by real-time quantitative reverse transcription-PCR (qRT-PCR) the expression of 17 circulating miRNAs, previously found to be significantly overexpressed in tissue pancreatic neoplasms, in a set of 182 plasma samples (94 PDAC, 19 IPMN, 18 chronic pancreatitis, and 51 disease-free controls). Then, we analyzed CA19.9 levels in the same plasma set, and we assessed the diagnostic values of differentially expressed miRNAs, CA19.9, and all possible combinations. RESULTS: Of note, 16, 14, and 9 miRNAs were significantly increased in PDAC, IPMN, and chronic pancreatitis, respectively, compared with control plasmas. miR-21-5p, miR-33a-3p, miR-320a, and miR-93-5p showed the highest discriminating capacity for pancreatic neoplasia (PDAC or IPMN) with an area under the receiver operating characteristic curve (AUC) of 0.86, 0.85, 0.85, and 0.80, respectively. 2-miRNA combinations improved these performances reaching AUC = 0.90 for "miR-33a-3p+miR-320a." Addition of CA19.9 increased the diagnostic potential of miRNA signatures even further achieving an AUC of 0.95 (93% sensitivity and 85% specificity) for the combination of "miR-33a-3p+miR-320a+CA19.9." CONCLUSIONS: Novel signatures combining miRNAs and CA19.9 could be used as noninvasive biomarkers for early detection of PDAC.

Integrated Analysis of Germline and Tumor DNA Identifies New Candidate Genes Involved in Familial Colorectal Cancer.

Colorectal cancer (CRC) shows aggregation in some families but no alterations in the known hereditary CRC genes. We aimed to identify new candidate genes which are potentially involved in germline predisposition to familial CRC. An integrated analysis of germline and tumor whole-exome sequencing data was performed in 18 unrelated CRC families. Deleterious single nucleotide variants (SNV), short insertions and deletions (indels), copy number variants (CNVs) and loss of heterozygosity (LOH) were assessed as candidates for first germline or second somatic hits. Candidate tumor suppressor genes were selected when alterations were detected in both germline and somatic DNA, fulfilling Knudson's two-hit hypothesis. Somatic mutational profiling and signature analysis were also performed. A series of germline-somatic variant pairs were detected. In all cases, the first hit was presented as a rare SNV/indel, whereas the second hit was either a different SNV (3 genes) or LOH affecting the same gene (141 genes). BRCA2, BLM, ERCC2, RECQL, REV3L and RIF1 were among the most promising candidate genes for germline CRC predisposition. The identification of new candidate genes involved in familial CRC could be achieved by our integrated analysis. Further functional studies and replication in additional cohorts are required to confirm the selected candidates.

High incidence of advanced colorectal neoplasia during endoscopic surveillance in serrated polyposis syndrome.

BACKGROUND: Serrated polyposis syndrome (SPS) has been associated with an increased risk of colorectal cancer (CRC). Accordingly, intensive surveillance with annual colonoscopy is advised. The aim of this multicenter study was to describe the risk of advanced lesions in SPS patients undergoing surveillance, and to identify risk factors that could guide the prevention strategy. METHODS: From March 2013 to April 2015, 296 patients who fulfilled criteria I and/or III for SPS were retrospectively recruited at 18 centers. We selected patients in whom successful clearing colonoscopy had been performed and who underwent subsequent endoscopic surveillance. Advanced neoplasia was defined as CRC, advanced adenoma, or advanced serrated lesion that were ≥ 10 mm and/or with dysplasia. Cumulative incidence of advanced neoplasia was calculated and independent predictors of advanced neoplasia development were identified. RESULTS: In 152 SPS patients a total of 315 surveillance colonoscopies were performed (median 2, range 1 - 7). The 3-year cumulative incidence of CRC and advanced neoplasia were 3.1 % (95 % confidence interval [CI] 0 - 6.9) and 42.0 % (95 %CI 32.4 - 51.7), respectively. Fulfilling both I + III criteria and the presence of advanced serrated lesions at baseline colonoscopy were independent predictors of advanced neoplasia development (odds ratio [OR] 1.85, 95 %CI 1.03 - 3.33, P  = 0.04 and OR 2.62, 95 %CI 1.18 - 5.81, P  = 0.02, respectively). During follow-up, nine patients (5.9 %) were referred for surgery for invasive CRC (n = 4, 2.6 %) or because of polyp burden (n = 5, 3.3 %). After total colectomy, 17.9 % patients developed advanced neoplasia in the retained rectum. CONCLUSIONS: Patients with SPS have a substantial risk of developing advanced neoplasia under endoscopic surveillance, whereas CRC incidence is low. Personalized endoscopic surveillance based on polyp burden and advanced serrated histology could help to optimize prevention in patients with SPS.

Mutational Signatures in Cancer (MuSiCa): a web application to implement mutational signatures analysis in cancer samples.

BACKGROUND: Mutational signatures have been proved as a valuable pattern in somatic genomics, mainly regarding cancer, with a potential application as a biomarker in clinical practice. Up to now, several bioinformatic packages to address this topic have been developed in different languages/platforms. MutationalPatterns has arisen as the most efficient tool for the comparison with the signatures currently reported in the Catalogue of Somatic Mutations in Cancer (COSMIC) database. However, the analysis of mutational signatures is nowadays restricted to a small community of bioinformatic experts. RESULTS: In this work we present Mutational Signatures in Cancer (MuSiCa), a new web tool based on MutationalPatterns and built using the Shiny framework in R language. By means of a simple interface suited to non-specialized researchers, it provides a comprehensive analysis of the somatic mutational status of the supplied cancer samples. It permits characterizing the profile and burden of mutations, as well as quantifying COSMIC-reported mutational signatures. It also allows classifying samples according to the above signature contributions. CONCLUSIONS: MuSiCa is a helpful web application to characterize mutational signatures in cancer samples. It is accessible online at http://bioinfo.ciberehd.org/GPtoCRC/en/tools.html and source code is freely available at https://github.com/marcos-diazg/musica .

Deciphering microRNA targets in pancreatic cancer using miRComb R package.

MiRNAs are small non-coding RNAs that post-transcriptionally regulate gene expression. They play important roles in cancer but little is known about the specific functions that each miRNA exerts in each type of cancer. More knowledge about their specific targets is needed to better understand the complexity of molecular networks taking part in cancer. In this study we report the miRNA-mRNA interactome occurring in pancreatic cancer by using a bioinformatic approach called miRComb, which combines tissue expression data with miRNA-target prediction databases (TargetScan, miRSVR and miRDB). MiRNome and transcriptome of 12 human pancreatic tissues (9 pancreatic ductal adenocarcinomas and 3 controls) were analyzed by next-generation sequencing and microarray, respectively. Analysis confirmed differential expression of both miRNAs and mRNAs in cancerous tissue versus control, and unveiled 17401 relevant miRNA-mRNA interactions likely to occur in pancreatic cancer. They were sorted according to the degree of negative correlation between miRNA and mRNA expression. Results highlighted the importance of miR-148a and miR-21 interactions among others. Two components of the Notch signaling pathway, ADAM17 and EP300, were confirmed as miR-148a targets in MiaPaca-2 pancreatic cancer cells overexpressing miR-148a. Moreover, a CRISPR-Cas9 cellular model was generated to knock-out the expression of miR-21 in PANC-1 cells. As expected, the expression of two miRComb miR-21 predicted targets, PDCD4 and BTG2, was significantly upregulated in these cells in comparison to control PANC-1.